Speed of sound in a superfluid Fermi gas in an optical lattice 
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A system of equal mixture of 6 Li atomic Fermi gas of two hyperfine states loaded into a cubic 
three-dimensional optical lattice is studied assuming a negative scattering length (BCS side of the 
Feshbach resonance). When the interaction is attractive, fermionic atoms can pair and form a 
superfluid. The dispersion of the phonon-like mode and the speed of sound in the long-wavelength 
limit are obtained by solving the Bethe-Salpeter equations for the collective modes of the attractive 
Hubbard Hamiltonian. 
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I. INTRODUCTION 

In the last decade the possibility of a superfluid alkali atom Fermi gas has attracted much attention both 
theoretically-^ - — and experimentally 1 '- because this phenomenon opens a new opportunity to study strongly cor- 
related quantum many-particle systems and to emulate high-temperature superconductors. Optical lattices are made 
with lasers, and therefore, the lattice geometry is easy to modify by changing the wavelength of the intersecting laser 
beams. Near the Feshbach resonance the atom-atom interaction can be manipulated in a controllable way because the 
scattering length a s can be changed from the BCS side (negative values) to the BEC side (positive values) reaching 
very large values close to resonance. We focus our attention on the BCS transition (negative scattering length) of 
degenerate fermionic gases to a superfluid state analogous to superconductivity. In particular, we consider an equal 
mixture of 6 Li atomic Fermi gas of two hyperfine states \F = 1/2, my = ±1/2 > with contact interaction loaded into 
an optical lattice. The two hyperfine states are described by pseudospins a =t,i- We also assume that the number 
of atoms in each hyperfine state per site (the filling factor) is smaller than unity, and that the lattice potential is 
sufficiently deep such that the tight-binding approximation is valid. The system in this case is well described by the 
single-band Hubbard model: 

H = -J 4>ia*Pj,<T - M n%,o + "i,1%4-- (!) 

Here, the Fermi operator C0i,a-) creates (destroys) a fermion on the lattice site i with pseudospin a =t,i and 
ni t a — ipj a ipi,a is the density operator on site i with a position vector r^. /i is the chemical potential, and the symbol 
2<y> means sum over nearest-neighbor sites. J is the tunneling strength of the atoms between nearest- neighbor 
sites, and U is the on-site interaction. On the BCS side the interaction parameter U is negative (the atomic interaction 
is attractive) . For simplicity we assume that each well of the periodic potential for atomic motion in three dimensions 
could be approximated by a harmonic potential. This harmonic approximation gives the following analytical results 
of J and U (H = 1)2: ' 

J = E R e ~ 

Here A is the laser wavelength, s is the lattice height, and m is the mass of the trapped 6 Li atoms. The recoil energy 
of the lattice En — ir 2 /2md 2 depends on the lattice constant d = A/2. In our numerical calculations the wavelength 
is chosen to be A = 1030 nm (Er = 1.293 x 10~ n eV)£. The lattice height is assumed to be s = 2.5. 

In what follows we study the spectrum of the collective modes of the Hamiltonian (JTJ . According to the Goldstone 
theorem, the long- wavelength limit of the spectrum has to be linear which means the speed of sound in this limit is 
independent of the wave-vector. In Ref. [6j the spectrum of the collective modes has been obtained from the poles 
of the density response function which had been calculated in the generalized random phase approximation (GRPA). 
This response-function version of the GRPA uses a 4 x 4 matrix L° (we follow the notations used in Ref. 6] ) which 
has nine (not six as it is stated in Ref. Jf|) independent elements: a, b, c, c, d, d, £1222, £2212 and £2222- Thus, 
the response-function version of the GRPA has produced incorrect expressions for the density response function (see 
Eqs. (26) and (27) in Ref. @). At zero temperature, the correct GRPA leads to the following Bethe-Salpeter (BS) 
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equations for the collective mode w(Q) and corresponding BS amplitudes G ± (k, Q) 16 : 



MQ)-e(k,Q)]G+(k,Q) = 

U bk,Q7q,Q + Kq 1 i,o\ G+(q, Q) - ^- Y [7k,q7q,Q - ^,q^,q] G"(q, Q) 
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Here the form factors are defined as follows: 7k, q = itk^k+Q + w k w k +Q, / k ,Q = UkUk+Q — I'ki'k+Q, 7k, q 



Mk^k+Q - Mk+Q«k, and m k ,Q = u k t; k+ Q + u k+Q t; k where u k = 1 



[1 + e(k)/S(k)] /2. The quantity e(k,Q) 



£J(k + Q) + E(k), where E(k) — ylS^+E? depends on the gap function A and the mean-field electron energy e k . 
We use a tight-binding form of the mean-field electron energy: e k = 2 J (cos k x d + cos k y d + cos k z d) — fi, where /i is 
the chemical potential. The gap function and the chemical potential have to be determined by the BCS number and 
gap equations: 
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where / = M/N is the filling factor, and we have M atoms distributed along N sites. 

The BS equations for the collective modes can be reduced to a set of four coupled linear homogeneous equations. 
The existence of a non-trivial solution requires that the secular determinant det|jx _1 — V\\ is equal to zero, where the 
bare mean-field-quasiparticle response function x & n d the interaction V = diag(—U, —U, U, —U) are 4x4 matrices: 
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Here we have introduced symbols I a ,b — F a: b(e(k, Q)) and J a ,b — F at b(ui), where F ai b(x) is defined as follows (the 
quantities a(k, Q) and b(k, Q) = l k ,Q, m k ,Q, 7k,Q or Yk.Cj): 



F a ,b{%) = T?Y/ 



xa(k, Q)6(k, Q) 



N ^ w 2 -e 2 (k,Q) 



It is worth mentioning that the GRPA equations for the collective mode derived by Belkhir and Randeria^i can be 
obtained if we neglect in ([5]) all elements with index 7. In this case x an d V are 3x3 matrices. 



II. SPEED OF SOUND IN A CUBIC LATTICE 



The velocity of sound is important because it tells us how fast the sound propagates in the system, but more 
importantly, it is intimately related to the normal (phonon) part of the liquid according to Landau's theory of 
superfluidity^. 

In our numerical calculations, the sum over k is replaced by a triple integral over the first Brillouin zone: 
— 7r < k x d < 7r, — 7r < kyd < 7r and — ir < k z d < 7r. After that, we applied the substitutions x — tan k x d/4,y — 

tan/cy<i/4 and z = tan/cjd/4 to rewrite the integrals in the form of Gaussian quadrature J dx dyf(x,y, z)/(l + 
x)(l + y)(l + z). The corresponding integrals are numerically evaluated using 49 x 49 x 49 (xi,yj,Zk) points: 
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FIG. 1: The dispersion w/A of the phonon-like collective mode. For filling factor / = 0.5, lattice height s — 2.5, and scattering 
length a s = — 1000as (as is the Bohr radius of hydrogen), the chemical potential /j, — 0.326.E.R and the gap energy A = 0.05-Er 
are obtained by solving the number and gap equations Q. The speed of sound in the long- wavelength limit is 8.1 mm/s. The 
puncture curve represents the dispersion calculated in Ref. |(|. 
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FIG. 2: The speed of sound (in units mm/s) as a function of the scattering length a s /d (d = 515 nm is the lattice constant). 
The filling factor is f=0.5, and the lattice height is s = 2.5. The puncture curve represents the speed of sound calculated in 
Ref. [fj. 

f_ t dx $_ x dy J_ 1 dzf(x,y,z)/(1 + x)(l + y)(l + z) = YntiYlj=iYlf=i w i w j w kf{xi,yj,Zk), where w l is the cor- 
responding weight. It can be checked that there is no difference between the approximation by integrals and the case 
when the sums over k are taken explicitly assuming 128 sites per dimension. 

In Fig. 1 and Fig. 2 we present the results of our calculations of the dispersion of the phonon-like mode and the 
speed of sound as a function of the scattering length assuming that the filling factor and the lattice height are / = 0.5 
and s = 2.5, respectively. The long-wavelength part of the dispersion is linear with sound velocity of about 8.1 mm/s. 
For higher momenta the dispersion saturates to 2A. As it is expected, when the interaction between the atoms is 
increased by increasing the scattering length, the compressibility of the system increases, and therefore, the speed of 
sound decreases, as can be seen in Fig. 2. In both figures, there exists a difference of about 10 -15 percents between 
the BS approach and the response-function calculations presented in Ref. [6]. 

III. CONCLUSION 

In this paper, we have used the BS equations in the GRPA to obtain the dispersion of the phonon-like collective 
mode and the corresponding sound velocity in the long- wavelength limit in the system of equal mixture of 6 Li atomic 
Fermi gas of two hyperfine states loaded into a cubic three-dimensional optical lattice. It is shown that the previous 
calculations, which have been obtained by studding the poles of the density response functions, are not in accordance 
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with our results derived by means of the BS equations in the GRPA. 
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